Top mass determination in leptonic final 
states with J/i/j 



X 

5-H 



Avto Kharchilava 

Institute of Physics, Georgian Academy of Sciences, Tbilisi, Georgia 



Abstract 

Q ■ Taking advantage of large top production rates at the LHC, the leptonic final states 

with J/tjj are explored for a precise determination of the top quark mass. The top 
is partially reconstructed by combining the isolated lepton from the W with the 
J/ip from the decays of the corresponding 6-quark. The method relies, to a large 
extent, on the proper Monte-Carlo description of top production and decay. The 
main emphasis is put on the expected systematics uncertainties. Mass measurement 
accuracy is dominated by the current understanding of theoretical uncertainties 
■ which result in a systematic error of ^ 1 GeV. 
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Ph. 1 Introduction 



One of the main interests in the top physics studies at the LHC will be an 
accurate measurement of the top mass. Being a fundamental parameter of the 
Standard Model (SM) it constraints, along with the W, the Higgs mass thus 
providing an excellent consistency check of the SM. 

Currently the best Tevatron single experiment results on the top mass are 
obtained with the lepton plus jets final states which yield: M top = 175.9 ± 4.8 
(stat.) ± 5.3 (syst.) (CDF) [1] and 173.3 ± 5.6 (stat.) ± 5.5 (syst.) (D0) [2]. 
The systematic errors in both measurements are largely dominated by the 
uncertainty on the jet energy scale and amount to 4.4 GeV and 4 GeV for the 
CDF and D0, respectively. On the other hand, the systematics in di-lepton 
channels are somewhat less, but the statistical errors are significantly larger, 
by a factor of > 2, as compared to the lepton plus jets final states. In future 
runs of the Tevatron with about 20-fold increase in statistics the top mass 
will be measured with an accuracy of ~ 2 GeV [3]; in the lepton plus jets 
channel the error is dominated by the systematics while in di-lepton channels 
the limiting factor would still be statistics. 
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The legitimate question to ask - what is the ultimate precision on the top 
mass that can be attained at the LHC? Apparently, to significantly improve 
on jet measurements at LHC would be a challenge. However, one can take the 
advantage of large top production rates and explore leptonic final states. In this 
paper we consider four- (three-) lepton final states involving the J/ip — > /x 4 fi~ 
decays, see Fig. 1. 



Fig. 1. Schematic of the top decay to leptonic final states with J/t/>. 

Here, the top mass is partially reconstructed by combining the isolated lepton, 
fi or e, from the W decays with the J/ip from the decays of the corresponding 
6-quark, i. e. reconstructing the lepton plus J /ij) invariant mass, M^j/^, [4,5]. 
In order to determine the top decay topology one can require an additional 
tagging muon-in-jet coming from other h. The branching suppression factor 
is 5.3 x 10~ 5 taking into account the charge conjugate reaction and W — > eu 
decays. In spite of rather low branching ratios, we stress that these final states 
are experimentally very clean and can be exploited even at highest LHC lu- 
minosities. Furthermore, one can also explore an other way to associated the 
J/-0 with the corresponding isolated lepton - by measuring the jet charge of 
identified fe('s) and not requiring the tagging muon. In this case we gain a 
factor of more than 10 in statistics as compared to the four-lepton final states. 

Similar approaches, namely the use of the correlation between the top mass 
and its decay products, e. g. in final states with an isolated lepton combined 
with the muon from the semileptonic decays of the corresponding fe's, have 
been explored in earlier studies [6,7]. One should say that all these methods 
of top mass determination essentially rely on the Monte-Carlo description of 
its production and decay. Nonetheless, the model, to a great extent, can be 
verified and tuned to the data. 



2 Analysis 

In the following we assume the ti production cross-section of 800 pb for M top 
= 175 GeV. Events are simulated with the PYTHIA 5.7 [8] or HERWIG 5.9 [9] 
event generators. Particle momenta are smeared according to parameteriza- 
tions obtained from detailed simulation of the CMS detector performance. 
Four-lepton events are selected by requiring an isolated lepton ofp T > 15 GeV 
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in central pseudorapidity range of \rj\ < 2.4 and three non-isolated, centrally 
produced muons of px > 4 GeV in \r]\ < 2.4, with the invariant mass of the two 
of them being consistent with the J/ip mass. These cuts would significantly re- 
duce the external (non-tt ) background, mainly the Wbb production, Q which 
can be further reduced by employing, in addition, two central jets from another 
W, if needed. Resulting kinematical acceptance of the selection criteria is 0.3; 
this rather small value is largely due to soft muons from J/-0 and b. In one 
year high luminosity running of LHC, corresponding to an L int = 10 5 pb _1 , 
and assuming the trigger plus reconstruction efficiency of 0.8, we expect about 
10 5 x 800 x 5.3 • 10~ 5 x 0.3 x 0.8 = 1000 events. 

An example of the £J/ip mass distribution with the expected background is 
shown in Fig. 2. The background is internal (from the ti production) and is 
due to the wrong assignment of the J/^ to the corresponding isolated lep- 
ton. These tagging muons of wrong sign are predominantly originating from 
B/B oscillations, b —* c — > \i transitions, W(— > c, r) — > \i decays, tt/K de- 
cays in flight and amount to ~ 30% of the signal combinations. The shape of 
the signal £J/i/j events (those with the correct sign of the tagging muon) is 
consistent with a Gaussian distribution over the entire mass interval up to its 
kinematical limit of ~ 175 GeV. The background shape is approximated by a 
polynomial function of degree 3. The parameters of this polynomial function 
are determined with "data" made of the wrong combinations oiiJ/i/) with an 
admixture of signal. In such a way the shape of the background is determined 
more precisely and in situ. Thus, when the signal distribution is fitted, only 
the background normalization factor is left as a free parameter along with the 
three parameters of a Gaussian. The result of the fit is shown in Fig. 2. We 
point out that this procedure allows to absorb also the remaining external 
background (if any) into the background fit function. 

As a measure of the top quark mass we use the mean value (position of the 
maximum of the distribution) of the Gaussian, M^y, . In four years running 
at LHC with high luminosity the typical errors on this variable, including 
the uncertainty on the background, are about 0.5 GeV. It is composed of 
< 0.5 GeV statistical error and < 0.15 GeV systematics contribution due to 
the uncertainty on the measurement of the background shape. 

The measurement of the M^j, can then be related to the generated top quark 
mass. An example of the correlation between the M max and M top is shown 
in Fig. 3 along with the parameters of a linear fit. For comparison, we also 
show the corresponding dependence in a more traditional isolated lepton plus 
//-in-jet channel. Not surprisingly, the stronger correlation, and thus a better 



No reliable event generator is currently available to simulate the Wbb process. 
PYTHIA results indicate that with the above cuts this source of the background 
can be kept at a per cent level. 
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Fig. 2. Example of the £J/ip invariant 
mass spectrum in 4-lepton final states. 
Number of events correspond to four 
years running at LHC high luminosity. 



Fig. 3. Correlation between M. max and 
the top-quark mass in isolated lepton 
plus 3/tp (solid line) and isolated lepton 
plus /x-in-jet (dashed line) final states. 



sensitivity to the top mass, is expected in the iJ/ip final states as compared 
to the isolated lepton plus //-in-jet channel. This is because, in the former 
case, we pickup a heavy object (i. e. J/tp ) from the 6-jet which carries a 
larger fraction of its momentum. Apparently, when determining the top mass, 
the M^]^ measurement error, statistical and systematic, would scale up as 
the inverse slope value of the fit which is a factor of 2 in our case. Hence the 
statistical error on the top mass in this particular example would be ~ 1 GeV. 

It is appropriate to comment on the ways to obtain a larger event sample. Pos- 
sibilities might include reconstruction of the J/ijj — > e + e~ decays in 6-jets [10] 
and/or relaxing the lepton thresholds [5]. An even larger event sample can 
be obtained in three lepton final states, i. e. without the tagging muon, but 
using instead the jet charge technique to determine the ti decay topology. The 
jet charge is defined as a pr-weighted charge of particles collected in a cone 
around the J/ijj direction. Obviously, this kind of analysis requires detailed 
simulations with full pattern recognition. However, particle level simulations 
performed with PYTHIA and with realistic assumptions on track reconstruc- 
tion efficiency indicate comparable to muon-tag performance figures, but with 
about 10 times less integrated luminosity. In any case, through the LHC life- 
time, one can collect enough events so that the overall top mass measurement 
accuracy would not be hampered by the lack of statistics; it would rather 
be limited by the systematics uncertainties which are tightly linked with the 
Monte-Carlo tools in use, as will be argued in the following section. 
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3 Systematics 

Let us first briefly mention the features of PYTHIA and HERWIG that are 
relevant to our analysis. They both incorporate only the leading order matrix 
elements for tt production and decay, and simulate the QCD final state effects 
via parton showers. This algorithm is not exact at higher orders in a s but 
could be, in principle, matched to the next-to-leading order calculations. Both 
event generators keep track of color connections, however the hadronization 
is performed in different ways: in PYTHIA (JETSET) it is done via a string 
fragmentation model while HERWIG exploits a cluster hadronization scheme. 
Eventually, free parameters of these event generators are tuned to reproduce 
the data, especially from e + e~ colliders. These tunings to data usually include 
event shape variables, identified particles production rates, their momentum 
spectra, etc. 

An essential aspect of the current analysis is to understand limitations which 
would arise from the Monte-Carlo description of the top production and de- 
cay. With PYTHIA we have investigated the following sources of systematic 
uncertainties on the measurements: 

(1) Initial State Radiation. 

The MJj a fy value is unchanged even switching the initial state radiation 
off. 

(2) Final State Radiation (FSR). 

Large shift of ~ 7 GeV is observed when it is switched off. To evaluate the 
uncertainty we varied the parton virtuality scale m m i n - an invariant mass 
cut-off below which the showering is terminated. A ±50% variation of it 
around the default (tuned to data) value of 1 GeV induces an uncertainty 
of t° i 5 GeV. 

(3) QED bremsstrahlung . 

Only a small effect is observed when it is switched off. Furthermore, QED 
radiation is well understood and can be properly simulated. 

(4) Parton distribution functions (PDF). 

Depending on the relative fraction of gluon/quarks versus x in various 
PDFs the top production kinematics might be different. No straightfor- 
ward procedure is available for the moment to evaluate uncertainties due 
to a particular choice of PDF. Here we compare results obtained with 
the default set CTEQ2L [11] and a more recent CTEQ4L [12] parameter- 
izations of PDFs. The observed change in the value is well within 

0.1 GeV. 

(5) Top pt spectrum. 

One does not expect significant uncertainties associated with the top 
Pt spectrum, i. e. its generated shape. This is primarily due to a good 
agreement between the parton shower Monte- Carlos and the next-to- 
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leading order QCD calculations over a wide range of pr up to 1 TeV. 
In Ref. [13] this comparison is done with the HERWIG event generator. 
Similar results we have obtained also with PYTHIA. However, to see an 
effect we have artificially altered the top px spectrum by applying a cut 
at the generator level. For example, a la effect shows up only when the 
Pt cut-off pushed up to 100 GeV. 

(6) Top and W widths. 

Kinematical cuts that one applies usually affect the observed Breit-Wigner 
shape (tails) of decaying particles. Conversely, poor knowledge of the 
widths may alter the generated iJ/ip mass spectrum depending on the 
cuts. Only little change in the MJ 1 ^ value is seen even with the zero- 
width approximation. 

(7) W polarization. 

A significant shift is found for the isotropic decays of W when compared 
to the SM expectation of its ~ 70% longitudinal polarization. In future 
runs of the Tevatron the W polarization will be measured with a ~ 2% 
accuracy [3], and at the LHC this would be further improved so that it 
should not introduce additional uncertainties in simulations. 

(8) ti spin correlations. 

A "cross-talk" between t and t decay products is possible due to experi- 
mental cuts. To examine this effect in details the 2 — > 6 matrix elements 
have been implemented in PYTHIA preserving the spin correlations [14]. 
No sizeable difference in the M^jyiJ value is seen compared to the default 
2^2 matrix elements. 

(9) b fragmentation, except FSR. 

As a default, in PYTHIA we have used the Peterson form for the b- 
quark fragmentation function with ej, = 0.005. Variation of this value by 
±10% [15] leads to an uncertainty of ^25 GeV. 

(10) Background. 

The uncertainty would be mainly due to an inaccurate measurement of 
the background shape and the systematics contribution of < 0.15 GeV 
quoted in previous section would scale down with increasing statistics. For 
example, already with ~ 10 4 events the induced uncertainty is < 0.1 GeV. 

(11) Detector resolution. 

Here we have considered only Gaussian smearing of particle momenta and 
the effect on the measurement uncertainty is negligible. A possible 

nonlinearity of the detector response can be well controlled having a huge 
sample of J/ip, Y and Z leptonic decays that will be available at LHC. 



A summary of these studies is given in Fig. 4. One sees an impressive stability 
of results for reasonable choices of parameters. The expected systematic error 
in the determination is < GeV which translates into a systematic 

error on the top mass of 5M top < j^f GeV. 
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Fig. 4. Summary on the systematics studies. 



Another way to estimate systematics uncertainties could be by a comparison 
of PYTHIA and HERWIG expectations. With HERWIG we have tried various 
tunings from LEP experiments as well as its default settings listed in Ta- 
ble 1. They all yield comparable results to each other and to PYTHIA results, 
and are within < 0.5 GeV. This corresponds to a systematic uncertainty 
6M top < 1 GeV. 

Table 1 

HERWIG parameter settings [9]. 



Parameter 


Default 


ALEPH II 


L3 


OPAL 


Aqcd 


0.18 


0.173 


0.177 


0.15 




0.75 


0.645 




0.7 


CLMAX 


3.35 


3.025 


3.006 


3.75 


CLPOW 


2. 




2.033 


1.3 


PSPLT 


1. 


0.984 




0.85 


CLSMR 


0. 




0.35 





2 This should be taken with a precaution - bugs have been found in the treatment 
of the top decays in HERWIG [16] - so its impact on current analysis should be 
evaluated once the improved version of HERWIG is publicly available. 
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One more feature specific to top decays is related to the matrix element cor- 
rections to parton shower simulations [16]. A significant shift is observed in 
the mean value of the isolated lepton plus B-hadron invariant mass spectra 
when these corrections are implemented in HERWIG [17]. On the other hand, 
this shift does not necessarily induces an additional systematic error (which 
is at least partially accounted for in 2 - the FSR effects) once it is prop- 
erly taken into account. Further studies are required when the corresponding 
implementation is publicly made available. 



4 Conclusions 



Abundant production of top quarks at LHC allows to use leptonic final states 
to precisely determine the top quark mass thus avoiding limitations that are 
ultimately associated to measurements with jets. 

• In £J/ip final states the top mass can be determined with a systematic 
uncertainty of < 1 GeV. 

• These final states are experimentally very clean and can be exploited even 
at highest LHC luminosities. 

• The method heavily relies on the Monte-Carlo description of top pro- 
duction and decay, however the model can be verified and tuned to the 
data. 

• The precision would be limited by the theoretical uncertainties which is 
basically reduced to the one associated with the t — > b transition. 

To conclude, this method of top mass determination looks very promising 
and particularly well suited for the CMS detector with its 4 Tesla field and 
ensuring precise lepton measurements. 
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